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Abstract 

Organic agriculture requires cultivars that can adapt to organic crop management systems without the use of synthetic 
pesticides as well as genotypes with improved nutritional value. The aim of this study encompassing 16 experiments was to 
compare 23 broccoli cultivars for the content of phytochemicals associated with health promotion grown under organic 
and conventional management in spring and fall plantings in two broccoli growing regions in the US (Oregon and Maine). 
The phytochemicals quantified included: glucosinolates (glucoraphanin, glucobrassicin, neoglucobrassin), tocopherols (8-, 
y-, a-tocopherol) and carotenoids (lutein, zeaxanthin, (3-carotene). For glucoraphanin (17.5%) and lutein (13%), genotype 
was the major source of total variation; for glucobrassicin, region (36%) and the interaction of location and season (27.5%); 
and for neoglucobrassicin, both genotype (36.8%) and its interactions (34.4%) with season were important. For 5- and y- 
tocopherols, season played the largest role in the total variation followed by location and genotype; for total carotenoids, 
genotype (8.41-13.03%) was the largest source of variation and its interactions with location and season. Overall, 
phytochemicals were not significantly influenced by management system. We observed that the cultivars with the highest 
concentrations of glucoraphanin had the lowest for glucobrassicin and neoglucobrassicin. The genotypes with high 
concentrations of glucobrassicin and neoglucobrassicin were the same cultivars and were early maturing Ft hybrids. 
Cultivars highest in tocopherols and carotenoids were open pollinated or early maturing F n hybrids. We identified distinct 
locations and seasons where phytochemical performance was higher for each compound. Correlations among horticulture 
traits and phytochemicals demonstrated that glucoraphanin was negatively correlated with the carotenoids and the 
carotenoids were correlated with one another. Little or no association between phytochemical concentration and date of 
cultivar release was observed, suggesting that modern breeding has not negatively influenced the level of tested 
compounds. We found no significant differences among cultivars from different seed companies. 
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Introduction 

Organic food consumption is in part driven by consumer 
perception that organic foods are more nutritious and simulta- 
neously less potentially harmful to human health [1-2]. Studies, 
such as Smith-Sprangler et al. [3], have concluded that there is 
littie evidence for differences in health benefits between organic 
and conventional products, but other studies have indicated that 
organic vegetables and fruits contain higher concentrations of 
certain plant phytochemicals associated with health promotion 
than those produced conventionally [4—8]. A number of these 
compounds are produced by plants in response to environmental 
stress or pathogen infection, providing a potential explanation of 
why concentrations of these compounds might be higher in plants 
grown in organic systems without application of pesticides [9] . In 
addition, higher phytochemical levels may be due to the effects 



that different fertilization practices have on plant metabolism. 
Synthetic fertilizers used in conventional agriculture are more 
readily available to plants than organic fertilizers [10]. Nutrients 
derived from organic fertilizers need to be mineralized, and the 
availability of these nutrients depends on soil moisture, temper- 
ature and level of activity of soil organisms [1 1] . Conventional 
systems seek to maximize yields, resulting in a relative decrease of 
plant phytochemicals and secondary metabolites [12-15]. Corre- 
spondingly, compounds such as phenolics, flavonoids, and indolyl 
glucosinolates may be induced by biotic or abiotic stress [16-17]. 

Broccoli is an abundant source of nutrients, including provita- 
min A (P-carotene), vitamin C (ascorbate), and vitamin E 
(tocopherol) [18]. It is also a source of phytochemicals associated 
with health benefits and these include glucosinolates, carotenoids, 
tocopherols, and flavonoids [19-21]. Verhoeven et al. [22], Keck 
and Finley [23] and Here and Biichler [24], reported that diets 
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rich in broccoli reduce cancer incidence in humans. Strong 
associations between consumption level and disease risk reduction 
exists for glucosinolates (anti-cancer), tocopherols (cardiovascular), 
and the carotenoids (eye-health) [25]. 

Sulfur containing glucosinolates are found in the tissues of many 
species of the Brassicaceae family. When glucosinolates are 
consumed, they are hydrolyzed into isothiocyanates (ITC) and 
other products that up-regulate genes associated with carcinogen 
detoxification and elimination. Aliphatic glucoraphanin (up to 
50% of total glucosinolates) and the indolylic glucosinolates, 
glucobrassicin and neoglucobrassicin are abundant in broccoli 
florets [20,19,26]. Glucoraphanin is hydrolyzed either by the 
endogenous plant enzyme myrosinase [27-28] or by gut microbes 
to produce sulforaphane, an ITC. The indole glucosinolates are 
tryptophan-derived in a similar but alternate biosynthetic pathway 
[29]. The health promoting effects of the indolyl glucosinolates are 
attributed to indole-3-carbinol, a hydrolysis product of gluco- 
brassicin, N-methoxyindole-3-carbinol and neoascorbigen, hydro- 
lysis products from neoglucobrassicin, and the catabolic products 
derived from alkyl glucosinolates. Clinical studies have shown that 
the glucosinolate hydrolysis products reduce the incidence of 
certain forms of cancer (e.g., prostate, intestinal, liver, lung, breast, 
bladder) [30-35]. The lipophilic phytonutrients found in broccoli 
include the carotenoids lutein, zeaxanthin, B-carotene, and 
tocopherols (forms of vitamin E) [36-37]. In addition to their 
role as vitamins, these compounds are powerful antioxidants [38- 
39]. Consumption of vegetables high in tocopherols and caroten- 
oids has decreased the incidence of certain forms of cancer [40] . 
Lutein and zeaxanthin protect against development of cataracts 
and age-related macular degeneration [41]. Tocopherols have also 
been associated with reduced risk of cardiovascular disease by 
preventing oxidative modification of low-density lipoproteins in 
blood vessels [42] . 

The genetic potential for high nutrient content has long been a 
concern of the organic industry in order to meet the expectations 
of organic consumers. This has often been manifested by 
questioning whether modern elite cultivars may have lower levels 
of nutritional content than older open pollinated cultivars. Indirect 
evidence supporting this argument comes from Davis et al. [43], 
who compared USDA nutrient content data for 43 garden crops 
released between 1950 and 1999. Statistically significant decreases 
were noted for six nutrients (protein, calcium, potassium, iron, 
riboflavin, and ascorbic acid), with declines ranging from 6% for 
protein to 38% for riboflavin. Crop varieties in 1950 had been 
bred to be adapted to specific regions and a relatively low input 
agriculture system, but contemporary cultivars are selected for 
yield, disease resistance, broad adaptation to high input agriculture 
systems, and for increased 'shipability' and shelf life. Traka et al. 
[44] recommend breeding with greater genetic diversity when the 
goal is enhanced phytochemical content by exploiting wild crop 
relatives. The genotype is important in determining the level of 
nutrients in a crop cultivar [45-47]. What is unclear, however, is 
whether the nutritional content of a cultivar is associated with 
certain genotypic categorization, e.g. old versus modern, open 
pollinated versus F[ hybrid cultivars. In addition, there is no clear 
differentiation as to what extent nutritional content in a crop is 
determined by genotypic or by field management factors or by the 
interaction of both. Some studies comparing performance of 
genotypes in organic and conventional production systems have 
shown that for certain agronomic traits, cultivars perform 
differently between the two production systems (e.g. for winter 
wheat: Murphy et al. [48], Baresel et al. [49]; for lentils: 
Vlachostergios et al. [50]; for maize: Goldstein et al.[51]), while 
others have shown no differences in ranking performance (for 



maize: Lorenzana and Bernardo [52]; for onions: Osman et al. 
[53]; for cereals: Przystalski et al. [54]). The results of these studies 
have profound implications for organic cultivar selection and 
breeding strategies and raise questions as to the need for cultivars 
to be bred with broad adaptability or specific adaptation for the 
requirements of regional organic production and for designing 
breeding programs that optimize phytochemicals in an adapted 
management system. 

Previous studies comparing organically versus conventionally 
grown broccoli for nutritional quality have been 'market basket' 
(off-the-shelf) studies [55-56]. Harker [57] explained that the 
limitation of market basket studies is that they either have 
purchased the products from the store shelf and cannot relate 
differences to specific growing conditions or that the number of 
cultivars is too small to generalize the results. While other studies 
have compared cultivars from one production season time period 
to another, knowledge of the actual cultivar and production system 
(soil quality, temperature, rainfall) was not available [58,43]. The 
concentrations and form of health-promoting nutrients in Brassica 
vegetables have been reported to vary significantly due to (1) 
genotype (cultivar and genotypic class) [59,20,26,60,21,37,61,44], 
(2) environmental conditions such as season [62-67], light [19], 
max/min temperature, irrigation [68-69], (3) genotype by 
environment interactions [19,70-71]; (4) management system 
including soil fertility [72-73], organic versus conventional 
[13,74-75], days to harvest [63-64], and (5) post-harvest 
management [76-77]. Identifying specific growing conditions 
and genotypes that produce cultivars with varying phytochemical 
content and putative disease-prevention activity could offer value- 
added commercial opportunities to the seed and food industry. 

In addition to research conducted on how broccoli genotypes, 
management system and environment interact for horticultural 
traits [78], we address in this paper the question of how do 
genotypes, management system and environment interact to 
determine the nutritional contributions of broccoli to the human 
diet. We studied the relative importance and interaction among 
genotypes (cultivars, genotypic classes) and environment {man- 
agement system [M: organic (O) or conventional (C)], season (S, a 
combination of year and season within year, i.e., fall 2006, spring 
2007, fall 2007, spring 2008), location (E)} in a set of 23 broccoli 
cultivars for floret glucosinolate, tocopherol and carotenoid 
concentrations grown under organic and conventional production 
systems in two contrasting broccoli production regions of the US: 
Oregon and Maine. Specifically we addressed the following 
questions: (1) what is the impact of organic management system 
compared to the environmental factors including climatic region, 
season and their interactions [Genotype (G) x Environment (E) x 
Management System (M)]? (2) is there a significant difference in 
phytochemical content between different genotypes and genotypic 
classes (old and modern cultivars; open pollinated and F! hybrid 
cultivars; early and late maturing cultivars; and between different 
commercial seed sources)? (3) what is the best selection environ- 
ment for a broccoli breeding program for enhanced phytochem- 
ical content? 

Materials & Methods 

Plant Material and Field Trial Locations 

Twenty-three broccoli cultivars including open pollinated (OP) 
cultivars, inbred lines, and F[ hybrids were included in field trials 
(Table 1). Cultivars were grown in a randomized complete block 
design with three replicates in Maine (ME)-Monmouth (Latitude 
44.2386"N, Longitude 70.0356°W); and Oregon (OR)-Corvallis 
(Latitude 44.5647"N, Longitude 123.2608"W)] with each location 
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including organically (O) and conventionally (C) managed 
treatments. Plots contained 36 plants, planted in three rows of 
1 2 plants at 46 cm equidistant spacing within and between rows. 
The 2006 trials had only 18 of the 23 entries, and the Oregon 
2006 trial had only two replicates at the organic location. Field 
trials were conducted for three consecutive years with one 
production cycle in Fall 2006, two production cycles in Spring 
and Fall 2007 and one production cycle in Spring 2008. The 
primary management differences between the organic and 
conventional field trial sites are outlined in Table SI in File 
SI, which describes the production system, soils, fertility 
applications, the applied supplemental irrigation, and weather 
conditions for the area of study. Further details of the field design 
are reported in Renaud et al. [78]. 

Field Data Collection 

As plots approached maturity they were evaluated three times a 
week for field quality and broccoli heads that had reached 
commercial market maturity (approximately 10 to 12 cm in 
diameter for most of the cultivars while retaining firmness). Field 
quality traits evaluated on a 1 to 9 ordinal scale included head 
color, bead size, and bead uniformity. Average head weight was 
determined by taking the mean of the five individual heads per 
plot. Head diameter averaged for five heads at harvest maturity 
from each plot. Maturity was based on days to harvest from 
transplanting date. Detailed procedures and horticulture trait 
performance data are reported in Renaud et al. [78] . 



Broccoli Floret Samples and glucosinolate, tocopherol, 
and carotenoid analysis 

In order to analyse nutritional compounds of the broccoli heads, 
the following procedure was followed: As plots approached 
maturity, five broccoli head tissue samples were harvested fresh 
from each subplot at each trial location and were composited into 
a single sample per replication. The samples were frozen at — 20°C 
and shipped in a frozen state to the University of Illinois, Urbana- 
Champaign where they were freeze-dried and assessed for 
nutritional phytochemicals. Each sample was analyzed for the 
glucosinolates (glucoraphanin, glucobrassicin and neoglucobrassi- 
cin), carotenoids (fS-carotene, lutein, and zeaxanthin), and 
tocopherols (8-, y-, a- tocopherol) by high-performance liquid 
chromatography (HPLC) analysis using analytical protocols 
described in Brown et al. [19] for glucosinolates, and Ibrahim 
and Juvik [37] for tocopherols and carotenoids. Glucosinolates in 
lyophilized floret tissue samples were extracted and analysed by 
HPLC using a reverse phase C18 column. Three hundred mg 
samples of broccoli floret tissue were weighed out for extraction 
and the HPLC quantification of the tocopherols and carotenoids. 

Statistical Analysis 

Various linear mixed models were used for the analysis of trait 
variation. We followed the same methodology as described in 
Renaud et al. [78], which was comparable to the approach 
followed by Lorenzana and Bernardo [52]. For fitting the linear 



Table 1. Overview of commercially available broccoli cultivars, showing origin, main characteristics, included in paired organic - 
conventional field trials 2006-2008. 



Cultivar 


Abbreviation 


Origin 


Cultivar Type* 


Date of Market Entry 


Maturity Classification* 1 


Arcadia 


ARC 


Sakata 


Pi 


1985 


L 


B1 10 


B11 


Rogers 


Pi 


1988 


M 


Batavia 


BAT 


Bejo 


Pi 


2001 


M 


Beaumont 


BEA 


Bejo 


Pi 


2003 


L 


Belstar 


BEL 


Bejo 


Pi 


1997 


L 


Diplomat 


DIP 


Sakata 


Pi 


2004 


L 


Early Green 


EGR 


Seeds of Change 


OP 


1985 


E 


Everest 


EVE 


Rogers 


Pi 


1988 


E 


Fiesta 


FIE 


Bejo 


Fi 


1992 


L 


Green Goliath 


GRG 


Burpee 


Pi 


1981 


M 


Green Magic 


GRM 


Sakata 


Pi 


2003 


M 


Gypsy 


GYP 


Sakata 


Pi 


2004 


M 


Imperial 


IMP 


Sakata 


Pi 


2005 


L 


Marathon 


MAR 


Sakata 


Fi 


1985 


L 


Maximo 


MAX 


Sakata 


Pi 


2004 


L 


Nutribud 


NUT 


Seeds of Change 


OP 


1990 


E 


OSU OP 


OSU 


Jim Myers, OSU 


OP 


2005 


E 


Packman 


PAC 


Petoseed 


Fi 


1983 


E 


Patriot 


PAT 


Sakata 


Pi 


1991 


M 


Patron 


PAN 


Sakata 


Pi 


2000 


M 


Premium Crop 


PRC 


Takii 


Pi 


1975 


E 


USVL 048 


U48 


Mark Farnham, USVL 


Inbred 


not released 


L 


USVL 093 


U93 


Mark Farnham, USVL 


Inbred 


not released 


M 



a Cultivar Type: F,: hybrid; OP: Open Pollinated; Inbred. 
b Maturity Classification: E: Early; M: Mid; L: Late. 
doi:1 0.1 371 /journal.pone.0095683.t001 
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mixed models, GenStat 15 (VSNi, 2012) was used. The models 
followed the set-up: 

y=E+R(E) + G + GxE+e. 

Here y is the phytochemical response. Term E represents the 
environment in a very general sense, it includes all main effects 
and interactions of Season (S), Location (L) and Management (M). 
For analyses per location, the terms involving L were dropped. 
Similarly, for analyses regarding a specific management regime, 
the terms involving M were dropped. Term R(E) is the effect of 
replicate within environment, and there were two or three 
replicates in individual trials. G and GxE axe genotype and 
genotype by environment interaction effects, respectively. Finally e 
is a residual. 

Variance components were reported as coefficients of variation, 
i.e., CV = WOjv/x with V the variance corresponding to specific 
effects and x the trait mean. Repeatability was calculated from the 
variance components in its most general form as H 2 = Vq/(Vq 
V GL /nL + + V GS /nS + V GM /nM + V GLS /(nL.nS) + V GLM /(nL.nM) 
+ V GSM /(nS.nM) + V GLSM /(nL.nS.nM) + VJ(nEnS.nM.nR)), 
where the variance components correspond to the terms in the 
mixed model above. The terms nL, nS, nM and nR stand for the 
number of locations (2: Maine and Oregon), number of 'seasons' 
(4: Fall 2006, Spring 2007, Fall 2007, Spring 2008), management 
(2; organic and conventional), and replicates (2 or 3). 

Genotypic means were calculated by taking genotypic main 
effects fixed instead of random in the mixed models above. 
Pairwise comparisons between genotypic means were performed 
using GenStat procedure VMCOMPARISON . Correlations on 
the basis of genotypic means were referred to as genetic 
correlations. Genotypic stabilities under organic and conventional 
conditions were calculated as the variance for individual genotypes 
across all trials in the system. 

To assess the feasibility of selection for organic conditions (the 
target environment) under conventional conditions, we calculated 
the ratio of correlated response (for organic conditions using 
conventional conditions), CR, to direct response (for organic 
conditions in organic conditions), DR, as the product of the 
genetic correlation between organic and conventional systems (r G ) 
and the ratio of the roots of conventional and organic 
repeatabilities (H c and H 0 respectively): CR/DR = r G H G /H 0 . 
A ratio smaller than 1 indicates that selection is better done 
directly under organic conditions when the aim is indeed to 
improve the performance in organic conditions. 

Results 

Comparison of phytochemicals means over the 
environments 

Glucosinolates. Across all trials, glucoraphanin levels were 
comparable between locations and seasons but were more variable 
at the individual location and season trial analysis level (Table 2). 
Glucoraphanin, glucobrassicin and neoglucobrassicin levels were 
comparable between organic and conventional treatments. Com- 
parisons of organic versus conventional by location and season for 
the glucosinolate phytochemicals are presented in Figure SI. 
Comparable levels of glucosinolates were observed in the organic - 
conventional comparisons within locations and seasons. 

Tocopherols. Across trials compared regionally, Oregon had 
higher levels of all three tocopherols compared to Maine (Table 2, 
Figure S2). The tocopherols 8- and y- were higher in Fall 
compared to Spring, but not so for a-tocopherol (Figure S2). 



Organic and conventional levels for all tocopherol concentrations 
were in the same range and not significantly different. When the 
three tocopherols were analysed by organic versus conventional 
within location and season, there were no clear significant 
differences in management system across the season and location 
combinations (Table 2, Figure S2). 

Carotenoids. Overall, Oregon had higher levels of lutein and 
P-carotene compared to Maine (Table 2, Figure S3) and 
comparative levels of zeaxanthin (Table 2, Figure S3). Spring 
produced higher levels of all carotenoids compared to Fall levels in 
contrast to the glucosinolates and the 8- & y- tocopherol 
concentrations. There were no significant differences between 
organic and conventional for any carotenoid measured. When 
carotenoids were analysed by management system within location 
and season, P-carotene showed significantly lower levels in Maine 
in the Fall compared to other location and season combinations 
(Figure S3). 

Partitioning of variance components 

Glucosinolates. For glucoraphanin across all trials in both 
regions, Genotype (G) main effect accounted for the largest 
proportion of variance, followed by GxLxS interaction (Table 3). 
There was no Management (M) main effect, but M contributed to 
the three (LxSxM and GxSxM) and four-way interactions 
(GxLxSxM). In contrast to glucoraphanin, Location (L) had the 
largest effect for glucobrassicin and neoglucobrassicin across all 
trials in both regions, followed by the LxS interactions. For 
neoglucobrassicin the S and G main effect was more important 
than for glucobrassicin. When trials were further partitioned by 
location, a G and S main effect was apparent for neoglucobrassicin 
in both locations; for glucobrassicin the S main effects was only 
apparent in Oregon and not in Maine (Table S2 in File SI). 
There was M main effect for glucobrassicin and neoglucobrassicin, 
but not for glucoraphanin, and no GxM interaction for all 
glucosinolates. 

Tocopherols. For 8- and y-tocopherol across all trials in both 
regions, the Season (S) main effect accounted for the largest 
proportion of variance (Table 3). In contrast the proportion of the 
variation associated with S for a-tocopherol across all trials was 
minor. For all three tocopherols there was minor to no M effect, 
but a large L main effect, being the greatest for y-tocopherol. The 
G main effect showed a similar pattern to L. 

Carotenoids. For all three carotenoids across all trials in both 
regions, the G main effect described a significant component of 
total variance and was of largest influence for lutein (Table 3). 
The S main effect played an important role for zeaxanthin, and to 
a lesser extent for lutein but not for P-carotene. For all three 
carotenoids the L effect was minor, but the LxS interaction for P- 
carotene was relatively large and mostly associated with Maine 
(Table S2 in File SI). There was no M main effect; only for P- 
carotene was there a small effect of the G xM interaction (mainly 
driven by Maine). 

Repeatability, genetic correlation and ratio of correlated 
response to direct response 

Organic versus conventional. In the present study, we were 
able to estimate the proportion of the genotypic variance relative 
to phenotypic variance, but because we did not have a genetically 
structured breeding population, we apply the term repeatability 
rather than broad sense heritability. Of the phytochemicals 
studied, repeatabilities for concentrations of seven of the nine 
were comparable or higher in organic compared to conventional 
systems (Table 4). Only for glucobrassicin and 8-tocopherol was 
repeatability under organic conditions lower than under conven- 
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tional. In the analyses 5- and oc-tocopherol had relatively low 
repeatabilities. The highest repeatabilities were for glucoraphanin 
(0.82-0.84), neoglucobrassicin (0.75-0.76), y-tocopherol (0.72- 
0.75), lutein (0.83-0.85) and zeaxanthin (0.76-0.77). Genetic 
correlations were high between organic and conventional for the 
glucosinolates, y-tocopherol and lutein (0.84—0.95), while 5- 
tocopherol, a-tocopherol, zeaxanthin and (3-carotene were lower 
(0.63-0.77). The ratio of the correlated response to direct response 
for selection in the organic system was less than 1 .0 for all traits. 

By location and season. For the glucosinolates, glucorapha- 
nin and glucobrassicin repeatability at each location, season and 
treatment trial were comparable and generally high (0.83-0.97) 
between organic and conventional trials, while no clear trend for 
neoglucobrassicin repeatabilities was observed between organic 
and conventional aside from being much lower than glucorapha- 
nin and glucobrassicin (Table S3 in File SI). For y- and ot- 
tocopherol, repeatabilities were comparable between organic and 
conventional, while for 8-tocopherol repeatabilities were compa- 
rable between systems or higher in conventional except for one 
paired trial. For the carotenoids, repeatabilities were comparable 
or higher in organic for all paired trials, while for lutein in seven of 
the eight paired trials organic was comparable or greater than 
conventional. Repeatabilities for zeaxanthin concentrations were 
comparable for six of the eight paired trials. 

Comparison of cultivar ranking for phytochemical 
concentration and stability across trials 

To determine trends in cultivars with both the highest 
concentration of phytochemical groups most stable across 
locations, seasons and production systems, phytochemical concen- 
trations were plotted against stability per genotype across trials. A 
group of cultivars were identified as both highest in concentration 
and most stable and are indicated in the highlighted 'red circle' per 
phytochemical (Figure 1A-I). For glucoraphanin, the same group 
of cultivars had both the highest concentrations and were the most 
stable across production systems (Figure 1A; Table S4 in File 
SI). While for glucobrassicin, a different set of cultivars had the 
highest concentrations across production systems (Figure IB; 
Table S5 in File SI). Overall stability of all cultivars across 
production system was less related to cultivar mean concentrations 
for glucobrassicin than for glucoraphanin. None of the cultivars 
with the highest concentration for neoglucobrassicin were in the 
top quartile for stability across trials; all cultivars with the highest 
neoglucobrassicin content were in the bottom half for stability 
(Figure IC; Table S6 in File SI). Some but not all cultivars that 
had the highest concentrations of a-tocopherol were among the 
top group for 8- and/ or y-tocopherol. There was no relationship 
between 5-tocopherol concentrations and stability, but both y- and 
Qt- tocopherols had higher concentrations associated with greater 
stability (Figure 1D-F; Tables S7-9 in File SI). Open 
pollinated and early maturing cultivars had the highest and most 
stable concentrations for all carotenoids. (Figure 1G-T; Tables 
S10-12 in File SI). 

Comparison of phytochemical concentration by 
genotype classification 

The open pollinated and F! hybrid cultivars were compared 
across trials for each phytochemical analysed (Figure 2A). The 
levels of glucoraphanin in Fj hybrids tended to be higher than the 
open pollinated cultivars. But the inverse trend was observed for 
glucobrassicin, which was supported by the ranking and stability 
analysis where the Fi hybrids showed higher levels and more 
stability across trials than the open pollinated cultivars for 
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Table 4. Repeatabilities, genetic correlation and ratio of correlated response to direct response for broccoli phytochemicals 
comparing organic versus conventional management systems over all trial season/location combinations, 2006-2008. 



Repeatability (H) 




c 


o 




C R Q rg / ^org 


Glucoraphanin 


0.84 


0.82 


0.84 


0.83 


Glucobrassicin 


0.70 


0.64 


0.88 


0.84 


Neoglucobrassicin 


0.75 


0.76 


0.94 


0.94 


8-Tocopherol 


0.50 


0.42 


0.73 


0.66 


y-Tocopherol 


0.75 


0.72 


0.95 


0.93 


a-Tocopherol 


0.23 


0.35 


0.61 


0.76 


Lutein 


0.83 


0.85 


0.93 


0.94 


Zeaxanthin 


0.76 


0.77 


0.77 


0.78 


P-Carotene 


0.62 


0.72 


0.63 


0.68 



a Average genetic correlation between conventional and organic production systems across locations. 
b Ratio of correlated response to direct response. 
doi:1 0.1 371 /joumal.pone.0095683.t004 
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Figure 1 . Broccoli cultivar stability from trials conducted in two locations over four seasons with two management systems plotted 
against phytochemical content. A. Glucoraphanin, B. Glucobrassicin, C. Neoglucobrassicin, D. S-tocopherol, E. y -tocopherol, F. ot-tocopherol, G. 
Lutein, H. Zeaxanthin, I. (3-carotene. See Table 1 for cultivar name abbreviations. The C or O at the end of the cultivar abbreviation indicates 
conventional or organic management system, respectively. 
doi:10.1371/journal.pone.0095683.g001 
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Figure 2. Mean phytochemical content of broccoli genotypic classes. A. Mean phytochemical content of broccoli F, hybrids versus open 
pollinated cultivars, and B. Mean phytochemical content of early, mid- and late-maturing cultivars grown across all trials at two locations (Maine and 
Oregon), in two seasons (Fall and Spring) and in two management systems (Conventional and Organic) and conventional management systems. See 
Table 1 for key to cultivar F1 hybrid versus open pollinated classification and maturity classification. Significance (* =P<0.05, ** =P<0.01,*** =P< 
0.001). 

doi:1 0.1 371 /journal.pone.0095683.g002 



glucoraphanin. The reverse was observed for glucobrassicin. For 
the carotenoids, the open pollinated cultivars had a significantly 
higher mean value of lutein and zeaxanthin and tended to be 
higher for (3-carotene compared to the Fj hybrids. 

Based on the results of our field trials, the 23 cultivars of 
broccoli were grouped into three distinct maturity classes: Early 
(55-63 days); Mid (64-71 days); and Late (72-80 days) and 
analysed for the effect of the maturity class on phytochemical 
content (Figure 2B). For glucoraphanin, late maturing cultivars 
had significantly higher content levels, while for the carotenoids, 
early maturing cultivars tended to have higher concentrations and 
were significantly higher for lutein. 

When cultivar performance between genetic material originat- 
ing from two primary broccoli breeding companies was compared 
for phytochemical content there were no significant differences 
with the exception of lutein, where company l's cultivars had 
significantly higher concentrations than those of company 2 (data 
not shown). 

A negative correlation between the date of release and levels of 
glucobrassicin (R 2 = 0.21; p = 0.03) (Figure 3) was observed, but 
no significant correlations for any other phytochemical were seen 
when 21 cultivars (the total set minus the two inbred lines) were 
analysed by their date of commercial release (1975-2005). 



Correlation analysis among phytochemicals and 
horticulture traits 

Phytochemical correlation across trials. Correlation 
among phytochemicals indicated that glucoraphanin was signifi- 
cantly negatively correlated to glucobrassicin (Table 5). Correla- 
tions between the glucosinolates and the tocopherols were not 
significant. Glucoraphanin and neoglucobrassicin were negatively 
correlated to all carotenoids but only lutein and glucoraphanin 
were statistically significant. Glucobrassicin demonstrated a 
positive trend with all carotenoids. No statistically significant 
correlations were observed within tocopherols. A-tocopherol was 
positively correlated, while y-tocopherol was negatively correlated 
to all carotenoids. There were no significant correlations for ot- 
tocopherol with carotenoids. All carotenoids were highly positively 
correlated with one another. 

Phytochemical correlation to horticulture traits across 
trials. A correlation analysis was conducted for six horticulture 
traits, derived from the field study component of this research, 
Renaud et al. [78], and the nine phytochemicals across trials. The 
results indicated that greater head weight and head diameter were 
significantly positively correlated with glucoraphanin and nega- 
tively correlated with glucobrassicin, 8-tocopherol and the 
carotenoids. Increasing days to maturity was positively correlated 
with glucoraphanin, and negatively correlated to carotenoids. 
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Figure 3. Regression of broccoli floret glucobrassicin concentrations on date of cultivar release for 23 cultivars grown across all 
trials in two locations (Maine and Oregon), in two seasons (Fall and Spring), in two management systems (Conventional and 
Organic), 2006-2008. 

doi:1 0.1 371 /journal.pone.0095683.g003 



Head color was significantly correlated with 8-tocopherol and the 
carotenoids, but not with glucosinolates or y- and a-tocopherol. 
Bead size and bead uniformity were positively correlated with 
glucoraphanin, neoglucobrassicin and y-tocopherol and negatively 
correlated with glucobrassicin and the carotenoids. 

Principal component biplot analysis: correlation between 
phytochemicals and cultivars by production system 

In the principal component analysis the first PC axis accounted 
for similar amounts of the total variation in both conventional and 
organic production systems (43.5% vs. 39.6%). The second PC 
axis showed a similar trend with 17.02% for conventional and 
16.93% for organic (Figures 4A and 4B). The first two PC axes 
together accounted for 60.53% and 56.57% of total variation for 
conventional and organic, respectively. The PCA biplot analysis 
supported our findings that carotenoids were highly associated 
across systems, while tocopherols were highly associated in 
conventional, but not in organic (tocopherols demonstrated the 
largest shift between production systems). Glucoraphanin and 
neoglucobrassicin were associated with one another, but not with 
glucobrassicin across production systems. Glucoraphanin was 
associated with a-tocopherol in organic, but not in conventional 
treatments. Glucobrassicin was associated with 8- and a-tocoph- 
erol in conventional, but not in organic treatments. 5-tocopherol 
had a higher association with the carotenoids in organic than 
conventional. The biplots show response of both cultivars and 
phytochemical traits to environment. Those cultivars close to the 
origin reveal little about the relationship of cultivars and trait 
vectors, whereas those located near the extremes of trait vectors 
are those with the highest (or lowest) values for those traits. 

Discussion 

Impact of organic management system compared to 
environmental factors on phytochemical content 

Few studies have specifically compared the levels of health 
promoting compounds in Brassica vegetable species grown under 
organic and conventional production systems [13,74-75]. To our 
knowledge, this investigation is the most comprehensive study with 



the broadest range of phytochemical compounds (9) and a diverse 
set of broccoli cultivars (23) over regions (2), and management 
systems (2), with Fall and Spring season trials (2 each). In this study 
organic versus conventional management systems contributed the 
smallest source of variation compared to genotype, region and 
season. Within the phytochemicals studied individual compound 
concentrations responded differently. All compounds showed 
genetic variation, but also a substantial proportion of variance 
components were accounted for by high level interactions 
(Table 3; Table S2 in File SI). While M main effect was 
generally small, it had a substantial contribution in three- and 
four-way interactions. In particular, many GxLxSxM interac- 
tions were large relative to other variance components. This 
indicates that for the phytochemicals, M did have an influence on 
G, but that there were no consistent patterns across locations and 
seasons that would have shown up as significant GxM. Rather in 
each season and location, the paired organic and conventional 
environments differed significantly from one another but each 
situation was unique. In contrast to many comparisons between 
organic and conventional production systems [79], it should be 
noted that in our trials, yields averaged over the years did not 
differ significandy between the organic and conventional man- 
agement systems [78]. 

Among the nine compounds, glucoraphanin was the most 
strongly influenced by genotype followed by lutein: supporting the 
findings of several other broccoli studies where variation in 
concentrations for glucoraphanin [19,70,65-66] and lutein [21,37] 
was primarily due to genotype. For y-tocopherol, genotype was a 
large source of variation, but this compound was equally 
influenced by location and season (also found by Ibrahim and 
Juvik [37]). For glucobrassicin and neoglucobrassicin the location 
was the largest source of variation, but also LxS interaction was 
very influential, particularly for neoglucobrassicin, which is 
supported by Kushad et al. [20] and Schonhof et at. [26]. 
Jasmonic acid, a signal transduction compound in plants, is up- 
regulated under conditions of plant stress, wounding, and 
herbivory. Increased endogenous levels or exogenous application 
of this compound (or methyl jasmonate) increases biosynthesis and 
transport of neoglucobrassicin to broccoli florets. This up- 
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A Principal components biplot (60.53%) B Principal components biplot (56.57%) 




GTO ATO GLB DTO 



PC-1 (43.51%) PC-1 (39.64%) 

Figure 4. Principal components biplot of phytochemicals (vectors) and 23 cultivars (circles) grown in four seasons in Oregon and 
Maine. A. Biplot for conventional production, B. Biplot for organic production. See Table 1 for cultivar name abbreviations. Trait abbreviations: GLR: 
Glucoraphanin; GLB: Glucobrassicin; NGB: Neoglucobrassicin; DTO: 8-tocopherol; GTO: y-tocopherol; ATO: ct-tocopherol; LUT: Lutein; ZEA: Zeaxanthin; 
BCA: (3-Carotene. 

doi:1 0.1 371 /journal.pone.0095683.g004 



regulation was not observed for glucobrassicin biosynthesis [17]. 
This could explain why neoglucobrassicin was primarily under the 
control of Location and L x S interaction in our study. Season was 
the largest variance component for 8-tocopherol and zeaxanthin, 
which contrasts with the work of Ibrahim and Juvik [37] who 
found genotype had the largest influence on these compounds, 
followed by genotype by environment interaction although this 
study was constrained by the fact that the experiment was 
conducted in only one location over two growing seasons. For the 
other compounds such as Qt-tocopherol and P-carotene, LxS and 
the GxLxS interactions were most important. 

Overall we found high genetic correlations between glucosino- 
lates in organic and conventional trials. When trial locations were 
analysed separately, M main effect was present for glucobrassicin 
and neoglucobrassicin. The mean concentrations of glucobrassicin 
and neoglucobrassicin in broccoli from Oregon organic trials had 
higher concentrations compared to Oregon conventional trials, 
while Maine trials were comparable between management systems 
(Table 2, Figure SI). These results can be explained by the 
larger environment effect on glucobrassicin and genotype by 
environment effect on neoglucobrassicin found in the variance 
component analysis indicating sensitivity of these compounds to 
abiotic and/or biotic stresses. Our location specific findings are 
supported by those of Meyer and Adam [13] who performed a 
comparative study of the glucosinolate content of store bought 
organic and conventional broccoli and determined that the indolyl 
glucosinolates, glucobrassicin and neoglucobrassicin were signifi- 
cantly higher in the organically grown versus the conventionally 
grown. Evaluation of 10 broccoli genotypes over two years by 
Brown et al. [19] further supports our findings and those of Rosa 
and Rodrigues [62], Vallejo et al. [64], and Farnham et al. [70], 
that variation in concentration for glucoraphanin was primarily 
due to genetic variation, while differences in glucobrassicin was 
due to environmental variation (e.g. season, temperature) and 
genotype by environment interaction. The significandy higher 



levels of glucobrassicin in Oregon in the Fall harvested trials 
compared to Maine could be attributed to the higher maximum 
temperatures and GDD in Oregon compared to Maine (Table SI 
in File SI). 

Compared to glucosinolates, there is substantially less research 
on the genotype by environment interaction of tocopherol and 
carotenoid phytochemical groups in broccoli, and no specific 
studies exploring the influence of organic production system. In 
our study, minor management system effect at the overall trial 
analysis level was observed for the tocopherols and for carotenoids, 
there was management system effect only for lutein in Oregon 
Spring trials. Picchi et al. [75] also did not find differences in levels 
of carotenoids in cauliflower in organic versus conventional 
systems. In the tocopherols, there were no significant differences 
in location, but for 8- and y- tocopherol concentration levels were 
higher in the fall compared to the spring, while for a-tocopherol, 
concentration levels were higher in the spring compared to the fall. 
For the carotenoids, there were no significant location differences, 
however there was a seasonal trend that all carotenoids were 
higher in spring compared to fall. Ibrahim and Juvik [39] found 
significant environmental variation among 24 broccoli cultivars for 
carotenoids and tocopherols which they attributed to the stressful 
production environments. Factors explaining the genotype and 
genotype by environment interaction components of variation in 
the carotenoids and tocopherols could be clarified by the fact that 
environmental stimuli are both up- and down-regulating genes 
associated with carotenoid and tocopherol biosynthesis. There is 
evidence in the literature that there are coordinated responses of 
the carotenoid and tocopherol antioxidants in vivo. There was a 
reduction in rape seed (Brassica napus) tocopherol content in 
response to increased carotenoid levels due to over expression of 
the enzyme phytoene synthase [80]. This response could explain 
the negative correlation between y-tocopherol concentration and 
the carotenoids observed in our trials. 
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Differences in phytochemical content between different 
genotypes and genotypic classes 

The partitioning of variance indicated that genotype was an 
important source of variation for all glucosinolates. The cultivar 
ranking and rank correlation analysis demonstrated that there was 
a pattern in genotype content of glucosinolates where cultivars 
with the highest concentrations of glucoraphanin had the lowest 
levels for glucobrassicin (Figures SI). In our trials, the range in 
glucoraphanin concentrations across cultivars was (1.15— 
7.02 umol/g DW, Table S4 in File SI), while glucobrassicin 
was 1.46-3.89 umol/g DW, Table S5 in File SI). Several of the 
cultivars with the highest concentrations of neoglucobrassicin were 
those that had the highest concentrations of glucobrassicin. Range 
in neoglucobrassicin concentrations across cultivars was 0.68- 
4.54 umol/g DW, Table S6 in File SI). In earlier studies, 
glucosinolate concentrations in broccoli have shown dramatic 
variation among different genotypes. Rosa and Rodriguez [62] 
studied total glucosinolate levels in eleven cultivars of broccoli and 
found ranges from 15.2-59.3 umol/g DW. Among 50 accessions 
of broccoli Kushad et al. [20] found glucoraphanin content ranges 
from 0.8-22 umol/ g DW with a mean concentration of 7. 1 umol/ 
g DW, while Wang et al. [61] found glucoraphanin content of five 
commercial hybrids and 143 parent materials ranging from 1.57- 
5.95 umol/g for the hybrids and 0.06-24.17 umol/g in inbred 
lines and Charron et al. [65] found ranges from 6.4-14.9 umol/g 
DW. While the means in our study are somewhat lower, they are 
within the range of other studies. 

A genotype effect was observed for tocopherols, but predom- 
inantly for y-tocopherol. The PCA biplots (Figure 4AB) and the 
correlation analysis (Table 5) demonstrated the high positive 
correlations between 5-tocopherol, a-tocopherol and the caroten- 
oids (a-tocopherol and P-carotene were also highly correlated in 
the Kushad et al. [20] study). The cultivar relationship to different 
phytochemicals was represented in the biplots as well as in the 
cultivar content and stability analysis (Figure 1). Many cultivars 
with the highest concentrations in the tocopherols and carotenoids 
were open pollinated cultivars, inbreds and early maturing, older 
F[ hybrids. Many of this same group were also relatively high in 
glucobrassicin concentrations. Kurilich et al. [38] found that 
carotenoid and tocopherol concentrations among 50 broccoli lines 
were highly variable and primarily genotype dependent. Specif- 
ically, levels of P-carotene ranged from 0.4-2.4 mg/100 g FW. 
Ibrahim and Juvik [37] also found broad ranges for total 
carotenoid and tocopherol concentrations among 24 genotypes 
ranging from 55-154 ug/g DW and 35-99 ug/g DW, respec- 
tively. Farnham and Kopsell [2 1] studied the carotenoid levels of 
nine double haploid lines of broccoli. Similar to our findings, lutein 
was the most abundant carotenoid in broccoli ranging from 65.3- 
139.6 ug/g DM. The sources of variation for lutein were 
predominantly genotype, followed by environment and GxE 
interaction, which also supports our findings. No genotypic 
differences were found for P-carotene in Farnham and Kopsell 
[21], which is in contrast to our findings. Overall, they found that 
most of the carotenoids measured were positively and highly 
correlated to one another as was observed in our study (Table 5). 
Kopsell et al. [81] found lutein levels in kale of 4.8-13.4 mg/100 g 
FW where the primary variance components for both lutein and 
P-carotene were also genotype and season. 

Our research aimed also to address the question whether the 
phytochemical content of broccoli cultivars is associated with 
certain genotypic classes, e.g. open pollinated vs. F! hybrids; older 
vs. newer cultivar releases; and between commercial sources. 
Broccoli is typically a cross-pollinated, self-incompatible crop 
species and cultivars are either open pollinated and composed of 



heterogeneous genetically segregating individuals, or F! hybrids 
produced by crossing of two homozygous inbred lines, resulting in 
homogeneous populations of heterozygous individuals. In the 
1960's virtually all broccoli grown was derived from OPs. By the 
1990's almost all commercial cultivars were hybrids [82]. 

In our trials with 18 F| hybrids (released between 1975-2005) 
and 3 open pollinated cultivars (released from 1985-2005), we 
found several interesting trends related to genotype and genotypic 
class performance as it related to the three groups of phytochem- 
icals. When analysing Fi hybrid and open pollinated cultivars, 
they also demonstrated different performance patterns depending 
upon the individual phytochemical or group of compounds 
analysed. When cultivars were ranked for content and stability 
per phytochemical, there were distinct trends for certain 
compounds such as late maturing, F! hybrids outperforming early 
maturing F] hybrids and open pollinated cultivars for glucor- 
aphanin, while the inverse was found for glucobrassicin and all 
carotenoids studied. This analysis was further supported by the 
PCA biplots that showed a strong relationship for select cultivars to 
certain phytochemicals or groups of phytochemicals such as 'OSU 
OP' to the carotenoids. When the full set of cultivars was divided 
into F! hybrid and open pollinated groups and the means 
compared by phytochemical, the results further supported the 
individual cultivar analysis where F! hybrids had higher mean 
values for glucoraphanin than the open pollinated cultivars 
(Figure 2A). Clear cultivar performance differences were 
identified where early maturing versus late maturing cultivars 
performed differently depending upon the phytochemical 
(Figure 2B). We also found that late maturing cultivars had 
higher concentrations for glucoraphanin than early maturing lines 
(and the inverse for glucobrassicin and the carotenoids). Picchi et 
al. [75] studied the quantity of glucosinolates of an early and late 
maturing cultivar of cauliflower grown in one conventional and 
three organic production systems, and found a significantly higher 
level of glucoraphanin in the later maturing cultivar compared to 
the early maturing cultivar in the organic production system. 
Another interesting trend was that cultivars with higher concen- 
tration levels for those phytochemicals whose expression is heavily 
influenced by environmental factors were not necessarily the most 
stable across trial environments; as was the case with neogluco- 
brassicin, 8- and y-tocopherol in our study. For traits where 
genotype played a more significant role in contributing to 
variation, cultivars with a higher concentration level tended to 
also be those that were most stable across environments as was 
seen for lutein and glucoraphanin concentrations. 

No significant differences were found for cultivar performance 
in phytochemical concentrations between genetic materials 
originating from two distinct commercial sources, with the 
exception of lutein (data not shown). When the full set of broccoli 
cultivars were analyzed for a correlation between date of release 
and mean level of phytochemical content across trials, no 
significant correlation was found with the exception of a negative 
trend for glucobrassicin (Figure 3). Our data does not support the 
idea that modern breeding for high yield performance and disease 
resistance necessarily leads to a trade-off in level of phytochem- 
icals. Previous reports examining the relationship between year of 
release and performance had focussed on wheat vitamin and 
mineral content [83-85], and mineral content in broccoli [86-87]. 
However these authors did not study phytochemical content and 
their results were equivocal on the question on an innate biological 
trade-off between increased yield and nutritional content. 

Not many studies have included two or more groups of 
phytochemicals. In our study with three phytochemical groups we 
found that phytochemicals demonstrating a negative correlation 
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with one another (e.g. glucoraphanin with the carotenoids), 
showed an inverse cultivar response: e.g. cultivars with highest 
concentrations of glucoraphanin were the lowest in the carotenoids 
and vice versa. When both horticultural traits and phytochemicals 
were analysed for their phenotypic correlation, head weight was 
significandy and positively correlated with glucoraphanin and 
negatively correlated with 8- and ot-tocopherol and the caroten- 
oids. Farnham and Kopsell [21] explained that negative correla- 
tions may occur as a result of increased biomass accumulation in a 
certain genotype that is not accompanied by increased carotenoid 
production, effectively lowering the carotenoid concentration in 
the immature broccoli florets when pigments are expressed. 
Comparatively, head color was highly correlated to the caroten- 
oids and negatively correlated to the glucosinolates overall. The 
cultivar 'OSU OP' was explicitiy bred for a dark green stem and 
head color, not only for a darker green dome surface but also for a 
dark green interior color between the florets of the dome and in 
the stem (personal communication, Jim Myers 2013). 'OSU OP' 
was the highest in overall carotenoid concentrations across trials as 
it is known that carotenoids are correlated with chlorophyll 
concentrations and the intensity of green pigmentation [88]. 

Perspectives on breeding broccoli for enhanced 
phytochemical content specifically for organic 
agriculture 

Our study included predominantly broccoli cultivars selected 
for broad adaptability in conventional production systems and not 
purposely bred for high phytochemical content nor for adaptation 
to organic agriculture. What we can conclude from our data is that 
there has been little change in levels of several phytochemicals 
over three decades of breeding. This may indicate genetic 
variation for phytochemicals is limited in elite germplasm, or it 
may be the result of the lack of selection tools for these traits. This 
may be changing with recent efforts to introgress high glucor- 
aphanin from B. villosa to produce the high-glucoraphanin Fj 
cultivar 'Beneforte' [89-90,44] . The seed industry needs to exploit 
known sources of variation in the genus Brassica to enhance levels 
of other health-promoting phytochemicals and to broaden the 
genetic diversity of commercial broccoli germplasm. Our finding 
of a strong correlation between dark green color and high 
carotenoid levels provides breeders with a simple and efficient 
means of increasing carotenoids. The three groups of phytochem- 
icals studied contribute to health promotion in different ways. As 
these groups are related to different metabolic pathways selecting 
for one compound does not necessarily inadvertently improve the 
other compounds, and may even result in negative correlation as 
we have seen in our data between glucoraphanin and the 
carotenoids. Although these compounds belong to different 
metabolic pathways, their production may be coordinated through 
regulatory feedback loops, or the structural and/ or regulatory 
genes controlling these pathways may be genetically linked. 

Designing a breeding program for broccoli high in glucosino- 
lates would require the following considerations generated from 
our research: (1) Glucoraphanin is a highly genetically determined 
compound with minor location and season main effects but with 
substantial GxLxS interaction. (2) Comparatively, glucobrassicin 
and neoglucobrassicin are more impacted by location and season 
and LxS interaction with highest glucobrassicin concentrations 
and largest range in our Oregon Fall trials and neoglucobrassicin 
highest in Oregon Spring trials. (3) Cultivar performance for 
glucoraphanin and glucobrassicin and neoglucobrassicin was 
negatively correlated indicating that there may be a trade-off 
between glucoraphanin on the one hand, and glucobrassicin and 
neoglucobrassicin on the other hand. (4) Selection for glucor- 



aphanin without consideration of horticultural traits would 
probably result in larger headed and later maturing cultivars. 
Conversely, selection for smaller headed, early maturing cultivars 
would favor glucobrassicin and neoglucobrassicin at the expense of 
glucoraphanin. 

A breeding program for broccoli for high tocopherol content 
would require: (1) Overall the tocopherols were more season, 
location and LxS dependent and had lower overall repeatabilities 
compared to the glucosinolates. In a structured genetic population 
where additive genetic variance could be partitioned, narrow sense 
heritability would likely be low, and increasing tocopherol content 
would best be conducted with breeding methods suited to low 
heritability traits; (2) 8- and y-tocopherols were both season 
dependent and fall grown broccoli had higher concentrations of 
these compounds across trials and a wider range of content levels, 
whereas levels of a-tocopherol were higher in spring but the range 
was comparable under both seasons. Thus, fall would be the 
preferred environment for breeding for these compounds; (3) there 
were no significant differences for location for 8- or y-tocopherol, 
but the average levels of ot-tocopherol levels were significandy 
higher in Oregon than Maine, suggesting greater potential for 
genetic gain in the Oregon environment. 

If the goal is to design a breeding program for broccoli 
enhancing the levels of carotenoids it would require the following 
considerations: (1) For all three carotenoids studied, genotypic 
variation, particularly for lutein, was relatively more important 
than location and season; (2) however, zeaxanthin exhibited a 
large S (spring) and LxS interaction. For both fS-carotene and 
lutein, spring grown broccoli had significandy higher levels than 
fall produced. Thus, selection for carotenoids would probably be 
more effective in spring than in fall; (3) early maturing and small 
headed cultivars had higher levels of carotenoids. Since most of the 
carotenoids are associated with the outer surfaces of the 
inflorescence, smaller broccoli heads with a greater surface area 
to volume ratio should show higher concentrations of these 
compounds; (4) because carotenoids have high G main effect good 
germplasm sources as indicated in Figure 1 have high concen- 
trations of carotenoids and demonstrated stability across environ- 
ments. As all three carotenoids are highly correlated with one 
another, selecting for one should effectively select for all; (5) 
selection for darker green colour more widely distributed 
throughout the tissues of the head should allow the breeder to 
relatively efficiently increase carotenoid content in broccoli. 

In closing, we want to address the question of selecting in an 
organic or a conventional environment. The argument commonly 
used to support selecting in productive environments is that 
heritabilities are higher compared to resource poor environments 
[91-92]. Organic is often considered a low-external input 
environment, resulting on average in 20% less yield compared 
to conventional production [79]. Nevertheless, in our trials 
repeatabilities for some phytochemicals were higher or compara- 
ble to conventional (Table 3). Narrow sense heritabilities would 
be expected to be significandy lower. For those traits where 
repeatabilities were higher or comparable, direct selection under 
organic systems could enhance selection gain. In all cases, the ratio 
of correlated response to direct response was less than one 
suggesting that direct selection would allow more rapid progress 
than correlated selection. Our data on phytochemicals did not 
show a wider range of levels under organic conditions as we found 
for horticultural traits in the same trials [78], however, in several 
cases, repeatabilities in organic production were higher than in 
conventional. 

To maximize efficiency in a breeding program, commercial 
breeders may seek to combine breeding for both conventional and 



PLOS ONE | www.plosone.org 



13 



July 2014 | Volume 9 | Issue 7 | e95683 



Broccoli Breeding for Organic Agriculture and Phytochemical Content 



organic markets, and a combination of strategies can be proposed. 
Some studies that utilized highly heritable (agronomic) traits, 
where cultivar yield performance ranked similarly between organic 
and conventional management systems and which had high 
genetic correlations, suggested that early breeding be conducted 
under conventional conditions, with the caveat that advanced 
breeding lines be tested under organic conditions for less heritable 
traits (e.g. Loschenberger et al. [93]; Lorenzano and Bernardo, 
2008) [52]. In studies where cultivar yield performance differed 
between management systems and there were significant differ- 
ences in cultivar ranking, and in some cases low genetic 
correlations for lower heritability traits (e.g. Kirk et al. [94]; 
Murphy et al. [48]), these studies recommended that cultivars 
intended for organic agriculture be selected only under organic 
conditions. In our study of phytochemicals, we would recommend 
for organic purposes selection under organic conditions for the 
compounds where genetic correlations between organic and 
conventional were moderate. 

Supporting Information 

Figure SI Comparison of broccoli cultivars for glucosi- 
nolates (umol/g DW) grown across all trials in two 
locations (Maine and Oregon), in two seasons (Fall and 
Spring), in two management systems (Conventional and 
Organic), and at the individual trial level, 2006-2008. A. 
Glucoraphanin, B. Glucobrassicin, C. Neoglucobrassicin. 
(TIF) 

Figure S2 Comparison of broccoli cultivars for tocoph- 
erols (umol/g DW) grown across all trials in two 
locations (Maine and Oregon), in two seasons (Fall and 
Spring), in two management systems (Conventional and 
Organic), and at the individual trial level, 2006-2008. A. 
8-tocopherol, B. y-tocopherol, C. ot-tocopherol. 
(TIF) 

Figure S3 Comparison of broccoli cultivars for carot- 
enoids (umol/g DW) grown across all trials in two 
locations (Maine and Oregon), in two seasons (Fall and 
Spring), in two management systems (Conventional and 
Organic), and at the individual trial level, 2006-2008. A. 
Lutein, B. Zeaxanthin, and C. P-carotene. 
(TIF) 

File SI Supporting tables. Table SI. Description of 
agronomic and environmental factors of the trial locations with 
paired organically and conventionally managed fields, 2006-2008. 
Table S2. Partitioning (%) of variance components for various 
traits of 23 broccoli cultivars grown across four pair combinations 
in Maine, season (Fall/Spring) and management system (Conven- 
tional/Organic), 2006-2008. Variance components reported as 
coefficients of variation. Table S3. Repeatability for broccoli for 
phytochemicals and per trial of 23 broccoli cultivars grown across 
eight pair combinations of location (Maine/Oregon), season (Fall/ 
Spring) and management system (Conventional/Organic), 2006- 
2008. Table S4. Glucoraphanin level (umol/g DW) of 23 
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